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We analyse the combined effect of the radiation-pressure and cross-Kerr nonlinearity on 
the stationary solution of the dynamics of a nanomechanical resonator interacting with an 
electromagnetic cavity. Within this setup, we show how the optical bistability picture in¬ 
duced by the radiation-pressure force is modified by the presence of the cross-Kerr interaction 
term. More specifically, we show how the optically bistable region, characterising the pure 
radiation-pressure case, is reduced by the presence of a cross-Kerr coupling term. At the 
same time, the upper unstable branch is extended by the presence of a moderate cross-Kerr 
term, while it is reduced for larger values of the cross-Kerr coupling. 
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I. INTRODUCTION 

In recent years, the exploration of the coupling between electromagnetic and mechanical degrees 
of freedom has witnessed unprecedented interest both from the applied and fundamental point of 
view (for a recent review see [1]). The physics of the optomechanical conplings has implications in 
fields as diverse as the investigation of gravitational wave physics through, e.g. the VIRGO |2] and 
LIGO [S] consortia, and the study of coherent quantum effects with a view to the manipulation of 
quantnm mechanical and optical degrees of freedom down to the level of single quanta. Prominent 
examples of the achievements connected to the manipulation of optomechanical degrees of freedom 
are represented by the gronnd-state cooling of mechanical resonators [1], the optical/microwave 
transdnction by mechanical means [5], the demonstration of potential quantum-limited amplifica¬ 
tion [6]. 

From this perspective, in cavity optomechanical systems strong emphasis has lately been put 
to the so called single-photon strong coupling regime, for which the dynamics of mechanical and 
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optical degrees of freedom are strongly coupled, leading e.g. to photon blockade effects [71 [8] or 
the appearance of multiple sidebands in the optical spectrum of the cavity [8] . 

Along these lines, it has recently been proposed [U] and experimentally demonstrated m that 
strong-coupling between single-photons and mechanical quanta of motion can be achieved in a 
microwave setup. In this setup, the addition of a strongly nonlinear inductive element (single- 
Cooper-pair transistor SCPT) allows for an increase by several orders of magnitude of the single¬ 
photon optomechanical coupling go. Beyond the pure enhancement of the radiation-pressure (RP) 
coupling strength, in [9] it has been shown how the inclusion of a nonlinear circuit QED element 
determines the appearance of higher-order terms in the effective interaction between optical and 
mechanical degrees of freedom. In this context, it has therefore become relevant to discuss how 
higher-order interaction terms modify the dynamical properties of optomechanical system (see 
e.g. m]), paralleling the analysis performed in different setups, such as membrane in the middle 
geometries |12j and, analogously, in experiments exploiting optomechanical coupling with ultracold 
gases daiii!. 

In this paper we discuss the role played by the lowest-order correction to the interaction term 
beyond the RP term: the so-called cross-Kerr (CK) interaction term. The CK term can be pictured 
as a change in the reflective index of the cavity depending on the number of phonons in the 
mechanical resonator, while the RP term depends on the displacement of the mechanical resonator. 
In recent years, the cross-Kerr interaction has been shown to play a crucial role in the field of 
quantum information processing; in particular in the definition of near deterministic CNOT gates 
m, for entanglement purification and concentration [ISlIITj, in the analysis of hyperentangled 
Bell states, and in the quantum teleportation of multiple degrees of freedom of a single photon 

HU. 

We focus here on how the CK coupling affects the stability of the system in presence of a strong 
driving coherent tone, focusing primarily on the optimally detuned condition on the red sideband. 
Our main goal is to show how the CK term plays a nontrivial role in the determination of the 
stability picture of the system under consideration, generally extending the parameters window 
within which the dynamics of the system is stable. 

II. THE MODEL 

In our analysis, the dynamics of the system is generated by the following Hamiltonian 


H — Ho + Hj-p + Hck, 


( 1 ) 
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FIG. 1. Schematic picture of the system, (a) Lumped element representation of the system in presence of 
SCPT (from [TU]). (b) Regular optomechanical setup, and (c) membrane-in-the-middle setup. 


where (hereafter h= 1) 


Hq = OJcCL^a + 


( 2 ) 


is the Hamiltonian for the non-interacting electromagnetic cavity/mechanical resonator system, 
where cOc and Wm denote the resonance frequencies of the cavity and the mechanics, and a (a^) and 
b (b^) are the boson annihilation (creation) operators for the cavity and the mechanics respectively. 
As it is usually the case in the analysis of optomechancial systems, we focus our attention on one 
cavity and one mechanical resonator modes. The RP coupling term Hj-p is given by 

Hrp = goa^a (b^ + b^ . (3) 

The term Hj-p can be understood as a shift of the cavity resonant frequency as a function of the 
displacement of the mechanical resonator x oc (6l-|-6). Depending on whether a Fabry-Perot cavity 
or a lumped element circuit picture is adopted, according to the standard cavity optomechanics 
description, the shift in the resonant frequency can be modeled either as a change in the physical 
length of the cavity or as a change in the capacitance induced by the displacement of the mechanical 
resonator (see Fig. [^. From both perspectives, the coupling constant go can thus be expressed as 

dojAx) 

Qo ^zpf 


dx 


( 4 ) 
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where Xzpt is the zero-point potion associated with the displacement of the mechancal resonator. 
In the lumped element description, the resonant frequency given by cOc = l/y/ LC{x), where, as 
mentioned above, we have allowed for the possibility of a mechanical-position dependence of the 
overall circuit capacitance, leading thus to the following expression 

coc dC{x) 

ffo = -x.pf — 

In [To], it has been proposed and experimentally verified that the introduction of a strongly 
nonlinear circuit element (the SCPT) will boost the optomechanical coupling by several orders of 
magnitude. Essentially the RP coupling between cavity and mechanics, in presence of the Josephson 
junction qubit, is mediated by the strong nonlinear inductance determined by the presence of the 
qubit. Furthermore, in [^ it was shown that the qubit-mediated coupling introduces higher-order 
terms in the coupling between optics and mechanics. In this paper we focus on the lowest-order 
coupling term relevant in the pump/probe setup, ubiquitous in the analysis of optomechanical 
systems; from |9|, invoking the rotating-wave approximation, this term is given by 

-f^ck = gcWab%, (6) 


where 


^ 


2 d'^UJc 


( 7 ) 


As previously mentioned, this term can be interpreted as an electromagnetic field phase shift 
induced by the number of mechanical phonons, as opposed to the RP term, which represents 
a phase shift of the field, depending on the displacement of the mechanical resonator from the 
equilibrium position. 

In our analysis, we describe the coupling between the optomechanical system and the environ¬ 
ment in terms of input-output formalism m- Within this framework, the equations of motion 
generated by the Hamiltonian Q for a and b can be written as 


a = —i 


+ bj + gckb'^b a - -a + 


Wc + go 

b = -i (^iVrn + gcka^a^ b - goo^a - 


( 8 ) 

(9) 


where k and 7 describe the coupling to the environment of the cavity and the mechanical degrees 
of freedom, respectively. 
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Adopting a standard approach, with a view to considering a situation in which the cavity 
is strongly driven by a coherent signal ain exp and in absence of a direct driving of the 

mechanical motion, we decompose the field operators o, 6, Om and bin so that a = a + 5a, b = /3 +6b, 
“in = ctin + b“in and bin = An + 6bin- Assuming that the amplitude of the fluctuations to be small 
with respect to the coherent signals, in a frame co-rotating with the input coherent drive at 
frequency Up, we can write the zeroth-order approximation to Eqs. Q and Q as 

/~Kain 


a = 


f - i [Ao - 50 iP* + P)- 5ck|^P] 


/3 = - 


5o|a| 


2 +i(a;ni + 5ck|aP)’ 


( 10 ) 

( 11 ) 


where we have introduced the detuning Aq = Wp — Wc and we have assumed An = 0. The solutions 
of Eqs. (10) and as a function of the input field amplitude ain, represent the equilibrium 
solutions for the cavity field amplitude a and for the amplitude of the mechanical oscillator /3, in a 
frame rotating at cOp. For 5ck = 0, Eqs. (10) and reproduce the steady-state solution obtained 
for a regular optomechanical system in presence of RP coupling mm- In that case, if the input 
field is detuned so that | A| = Wm allowing for the optimal exchange of energy between cavity and the 
mechanical modes, in particular Aq = —Wm represents the optimal cooling condition (see mEH), 
while Aq = Wm allows for optimal amplification of an incoming signal around the cavity resonant 
frequency [6|. It has recently been shown [22| that, in presence of CK coupling and for moderate 
driving, such as to allow for the use of pure RP steady-state solutions, the optical damping of the 
mechanics is a non-monotonous function of the optical drive, thus potentially hindering the cooling 
of the mechanical motion by optical means for a red-detuned pump, and on the other, limiting the 
parametric instability for a blue sideband drive. 

In the present work we discuss a different aspect of the problem, by considering the stability 
properties of the solutions of Eqs. (10) and ©• The determination of the stability of the system 
is performed through the standard stability analysis of linear autonomous systems of differential 


equations, allowing us to identify the stability character of the equilibrium solutions of Eqs. (10) 
and (0 as a function of the different parameters characterising the system, with particular focus 
on 5ck- 


III. STEADY-STATE SOLUTIONS 


From Eqs. (10) and it is clear how the CK interaction can modify the optical bistability 
picture associated with the RP only setup (Fig. [^: due to the presence of higher-order terms both 
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in a and (3, in principle, a more complex stability diagram as a function of the external drive should 
be expected, more specifically a larger number of (stable or unstable) equilibrium points should 
appear. However, while some important qualitative and quantitative differences do arise because 
of the presence of the CK term, no extra physical equlibrium solutions are actually present. More 


specifically, combining Eqs. (10) and (11), we obtain a quintic equation for |ap, which represents 


the mean-field cavity occupation - we note here that, for pure RP {gck = 0), the equation would be 
cubic. Moreover, both for the pure RP case and in presence of CK coupling, lap can have either 
one or three positive and real solutions, meaning that in the quintic equation associated with the 
CK case, two of the solutions violate the condition |ap > 0 for all values of the parameters, and 
thus have to be discarded. 


For positive detunings Aq > 0, all values of |ap corresponding to the solutions of Eqs. (10) 


and © would be complex or negative, therefore implying that the optical bistability can only be 
found for negative detunings (red sideband, see also [20]). The values of a^ and ac, along with 
the corresponding input-field values and points define the branch B — C in Fig. (a). As 
previously mentioned, in this case, the presence of the nonzero CK term, due to considerations 
concerning the unphysical nature of some of the solutions to eqs. (10) and ( [IT| ), does not introduce 
further equilibrium points in the dynamical system analysis. However, as it is possible to see from 
Fig-d the presence of a CK coupling term shifts the relative position of points B and C, reducing 
the range of pump amplitude aiN for which the branch B — C is present. In order to discuss 
the stability character of the steady-state solutions, we investigate the stability properties of the 
following dynamical system 


6a = iA6a — —6a + iG{6b + 6b^) -|- \^6ain 

6b = -ioj^6b - ^6b + i{G*6a + G6a^) + ^6bi^ (12) 


which represent the first-order expansion of Eqs. ([^ and ([^ around the equilibrium points defined 
by Eqs. and and are formally analogous to the first-order equations obtained for the 

pure RP case (see e.g. [T]). The essential difference between the RP and the CK case, is that the 
linearised coupling G, the effective detuning A and the effective mechanical frequency ui^ are given 
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by the following expressions 


G — goa ( 1 + 


9ck I Oi I 


Wm + £/ck|ap 


A = Ao + 2^|a 


9o 


Wn 


1 + 


5ck I cr I 


(Wm +5ck|aP)^- 


“1“ 9ck |c^| 


(13) 


IV. STABILITY ANALYSIS - LOWER BRANCH 


In this section we study the stability character of the solutions of the cavity optomechanical 
system, within the usual framework of stability analysis for autonomous linear systems. More 
explicitly, the nonlinear differential equations system given by Q and can be solved perturba- 
tively order-by-order. The solution of the zeroth-order equation corresponds to the steady state 
solution given by Eqs. (10) and while the first-order equation of the fluctuations around the 
stationary solutions, given by Eqs. (fl^, can be written as 


Sv = Av + 5vin 


1 T 


where Jv = [(5a, 56, (56'^] , 5 vin = 5aiN,5o|j.,f , 56 in, , and A is given by 


A = 


iA-l 


-iA - f -iG^ 


iG 
-iG* 
? 0 


2 


(14) 


(15) 


0 iG 

0 

iG iG* 

-iG* -iG 0 

The stability character of the steady-state solution is then provided by the sign of the solutions 
of the characteristic equation associated with A, pa(A) = 0. It is possible to show how, from the 
Routh-Hurwitz criterion, the stability of the system is characterised by the two following conditions, 

y^A^ 


+ 


^2 e 2 

^ + 4^2 A< + A 2<2 ^ ^ Q 


4 ' 16 ' - m ' m ■ 4 

[(7 -h K)ri - r2] r2 - (7 -h k)'^T3 > o 


(16) 

( 17 ) 


with 


Tl — ^+A2-|-7K-|- — -I- a;^2 

r 2 = ^+7A2 + ^ + K <2 


r3 = 


7 


2 a 2 


K-2,,,e 2 

+ + 4A|G|2^S. + 


4 

















8 


(a) 



M. 


aiN 


(b) 



FIG. 2. Optical bistability in the RP only case. The unstable branch B —C corresponds to the lower branch, 
associated with the optical bistability picture (condition ([I^). The unstable branch D — E corresponds to 
the violation of condition ([l7|. Plot parameters: Aq = —Wm, go = k = O.bwm, 7 = 0.12wni 


which correspond to the conditions oq > 0 and 030201 — {af + o|ao) > 0 , for the characteristic 
polynomial PA (A) = A'^ + 03 A^ + 02 A^ + oiA + oo, the other conditions for the stability of the system 
being identically satisfied for the system under consideration. The violation of the condition oq > 0 
corresponds to the branch B — C . As outlined in [20], in the pure RP case, up to linear order in 
y/cjm, the critical values for |ap, corresponding to the points B and C (see Fig. are given by 


I |2 ^0 

\OiB,c\ = -T. - 



1 - 


3k2 

4 ^ 


( 18 ) 
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indexes B and C correspond to the lower and upper sign respectively. In presence of CK coupling, 
the stability picture is modified. For small values of the presence of a CK coupling increases 
the range of values of |ap for which the branch B — C is unstable. In the resolved-sideband regime 
(k <C ujm = l^ol) the CK result is obtained as a (first-order) expansion in from the pure RP 
solution 


l«tl" = - 


^0 


|^ck|2 

l®C I 


Ao 




_ gckAp 

IGAoWm 

S'ckAo 


K 


9 3k2 \ 

(19) 

+ -t) 

9\ 

l> 

0 to 

1 

(20) 


2lJra ' IGAoWm 

However, increasing ^^k further, leads to a reduction and, eventually, the disappearance of the 
unstable branch B — C (see Figs. |^and[^. In the regime for which ^cklo^P 1) the critical value 
for which the unstable branch B — C disappears can be approximated by 


9ck = 


8ioi 


16- 


T 

4a;2 


( 21 ) 


V. STABILITY ANALYSIS - UPPER BRANCH 


The condition given by Eq. (16), derived from the condition 030201 —(af-|-a|ao) > 0, corresponds 
to the unstable branch D — E. For the pure RP case, the endpoints of the D — E branch are (up 
to linear order in 'j/k and Aq = —Wm) 

~ ( 22 ) 


I i2 —9 

I«d| = 5o 


I |2 -2 

= Go 



2e 4 32 16 


(23) 


(where e = 'j/k) Analogously to the approach employed to evaluate the effect of the CK term on 
the B — C branch in the small gck limit, we consider the perturbative correction to the RP result 
induced by gck- Up to linear order, the boundaries of the unstable region D — E, for small values 
of gck ( again for y/Zc <C 1) are given by 

{i + m? 


^ck 

Od 


^ ck 

Oe 


2g'^{l + SgB + vl) 
2 3(l + r/E)2 


325i^Ae(1 -I- 2r]Ey vV 


3k 


+ 3 


(24) 

(25) 


. _ 3(1 -b t?d,e)(1 + 2t/d^e) - ??d,e 

8(1-I-r/D,E)(l + 277 d,e)^ 

Qc]< I I o 

^D,E = —^k^D,E • 

% 


with 
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Within the validity limit of the perturbative approximation, the unstable region D — E is 
extended by the presence of the CK nonlinearity. Our numerical analysis shows that the region 
D — E increases with increasing until a critical value g*]^, for which E —>• oo, implying that, for 
^ck = 5 ck 5 upper branch is unstable for every |q;| > |q;d|- The value of can be evaluated 
considering, in the limit |ap 3> 1, the condition given by Eq. ( |17[ ) is expressed by a fourth-order 
polynomial in \a\^, Q{\a\^). Furthermore, the asymptotic behaviour of Q{\a\^) can be deduced 
from the sign of the highest-power coefficients € 3 ° and . In this limit, we have 


C 3 — 


32 


Wn 


, S'ck f 5ck7«^Ao 


OJt, 


cf=-iK 


9o 


- 2 

% 

2 


I 45 , 




+ (7 + k ) 


(26) 

(27) 


If we assume 2Aq > g^i^, it is possible to see how < 0 for every value of the cavity field, while 
c“ > 0, except when gc^ = 2g^l<jj^. In this case IoeI^ —t 00 , and the system remains unstable for 
every value of the cavity field above | od | • If ffck is further increased above g\ = 2gQ /Wm the region 
D — E becomes finite again and eventually disappears for g^. In the limit 7 /^ <C 1, 


2 2(7 -I- k)^ 2 r, * 

- - L/\ go ^ + 77 ^ 

7A^Ao tOml 


Aq ——OJn 


(28) 


VI. CONCLUSION 

In this paper, we have shown how the addition of an additional cross-Kerr coupling term to 
the usual radiation-pressure coupling term in the description of the dynamics of optomechanical 
systems substantially modifies the stability properties of the corresponding semiclassical equations 
of motion. In particular, we have shown how the inclusion of a CK term leads to an alteration and, 
eventually, to the disappearance of the unstable branch B — C, for ^fck > ffOk. Analogously, the 
unstable branch D — E is affected by the presence of a CK coupling term: for values of gs^ G ]0, 
the branch D — E is extended (E —)• 00 , for g'^k = 9c\i )> while, upon further increase of ^cki the 
branch D — E shrinks and, eventually, disappears for g^^ = g^^ 

cross-Kerr coupling constant gch is increased. 
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FIG. 3. Effect of gck on the stability diagram as a function of the input drive, weak ^ck- (a) The presence 
of the CK term reduces the input drive range for the branch B — C (increasing ^ck for the different curves 
from left to right), (b) The CK term first induces a broadening of the upper unstable region, which, upon 
further increase, is reduced (increasing gck for the different curves from bottom to top). Plots parameters, 
9ciilgo = 0,0.2,0.4; other parameters same as in Fig. 
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FIG. 4. Effect of ^ck on the stability diagram as a function of the input drive, strong g^k- (a) The presence 
of the CK term reduces the input drive range for the branch B — C (increasing ^ck for the different curves 
from left to right), (b) The CK term first induces a broadening of the upper unstable region, which, upon 
further increase, is reduced (increasing gck for the different curves from bottom to top). Plots parameters, 
gck/go = 0, 0.95,1 other parameters same as in Fig. 
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FIG. 5. Same as in Fig. for the puer RP case (gck = 0), and for ^ck = 0.2. The vertical dashed lines 
correspond to the input drive values for which the Ag-dependence is depicted in Figs. [6]|^ 
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(a) 


(b) 



Aq 





Aq 




FIG. 6. Ao dependence of the cavity population in (a) the pure RP case and (b) in the CK (^ck = 0.2) case. 
Increasing oin = 0.2,0.28,0.35 in units of Wm/ffo for increasing peak height. The presence of the CK term 
leads to a larger parameter range for which stability is observed. 
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FIG. 7. Aq dependence of the cavity population in the pure RP case. Increasing oin = 2,4,10 in units of 
w/go for increasing peak height (see inset). The case oin = lOw/go outlines the mechanism thorough which 
the result for Aq = —oj is stable at large drives (vertical dashed line). 
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FIG. 8 . Aq dependence of the cavity population in the CK (^ck = 0.2) case. Increasing aiN = (2,4, lOw/go) 
for increasing peak height. It is possible to see how, for large values of the input field, the system is stable 
for every value of detuning Aq. The case aiN = Scj/^q outlines the peculiar (“reentrant”) peak shape for 
intermediate input drive. 









